function rest = coefficientC(theta, beta, p)

rest = sqrt(binomial(theta,beta) ...
    *binomial(theta+p,p)*binomial(theta-beta+p,p)); 